Gamma line radiation from stellar winds in Orion complex: 
A test for cosmic ray origin theory 
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ABSTRACT 

We consider 7-ray emission from heavy nuclei that are accelerated in the shocks of 

stellar winds and are excited in interactions with the ambient matter in the Orion complex. 

We show that the recent detection of 7-ray lines in the Orion complex by COMPTEL 

telescope can be explained by such a scenario, assuming cosmic abundances. The scenario 

is consistent with recent models of cosmic ray acceleration in stellar winds of massive stars. 
O" 
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Bloemen et al.(1994) have recently reported the detection of 7-rays in 3 — 7 MeV range 



in the direction of the Orion complex with the COMPTEL telescope. The flux observed in 
that energy range is ~ 10 -4 photon cm -2 s . Subtracting the background radiation in a 
3° region covering Orion A and B, they suggested that the observed photons are dominated 
by lines at 4.4 and 6.1 MeV, and tentatively identified them with the de-excitation lines 
of 12 C* and 16 0*, respectively (also see, Cameron 1994). 

That these 7-ray lines should be emitted in interactions between high energy particles 
and interstellar nuclei, was anticipated long ago (see, for example, Meneguzzi and Reeves 
1975; Ramaty, Kozlovsky and Lingenfelter 1979). Bloemen et al. (1994) claimed that the 
observed line width is as large as ~ 1 MeV and inferred that the lines result from collisions 
of energetic C and O nuclei with ambient protons or a nuclei (rather than collisions of 
low energy cosmic ray protons or a-particles with ambient C and O nuclei in the gas and 



dust). They, however, doubted that standard abundances of C and O nuclei in cosmic rays 
could explain the observed luminosity. Bykov and Bloemen (1994) have recently invoked 
the overabundance of oxygen nuclei that is predicted in supernovae ejecta of massive stars 
to explain the observed 7-ray photons. 

In this Letter we show that cosmic ray C and O nuclei, accelerated in the shocks 
in stellar winds from massive O and B stars, can explain the observed 7-ray luminosity, 
assuming cosmic abundances. Shocks in the stellar winds of massive stars have been 
previously suggested as the accelerators of cosmic rays to moderate energies to energies as 
high as 3 x 10 9 GeV (Biermann 1993; Biermann and Cassinelli 1993; Stanev, Biermann, 
Gaisser 1993). Such shocks have also been claimed to generate X-rays (MacFarlane and 
Cassinelli 1989) and 7 rays from n° decays (White 1985; White and Chen 1992). Below 
we show that the COMPTEL observations can be understood in terms of these shocks in 
the stellar winds that abound in the OB associations of the Orion complex. 

2. STELLAR WINDS AND 7-RAYS 

The scenario for the production of 7-ray line photons is briefly sketched as follows: 
heavy nuclei are accelerated in the shocks in the stellar winds of OB stars and these nuclei 
are excited to emit 7-ray photons when they interact with cool gas surrounding these stars 
or the OB associations. In Orion A, for example, the OB association is embedded in a hot 
region which interacts with the surrounding dense and cool molecular gas whose density 
is ^j 10 5 cm -3 and kinetic temperature is T ~ 100 K (the Orion-KL region; Genzel and 
Stutzki 1989). Below, we will calculate the 7-ray luminosity expected in this scenario. 

The cross-section of the interaction between carbon nuclei and protons that leads to 
the emission of a 4.4 MeV photon by the excited carbon nuclei has been discussed by Ra- 
maty, Kozlovsky and Lingenfelter 1979. They further noted that the 4.4 MeV level in 12 C* 
can be more efficiently populated by spallation interaction, namely, cr[ 16 (9(p, 74.4) 12 C] ~ 
10~ 25 ( 25I g eV )~ - 6 cm" 2 for E £ 25 MeV (Ramaty, Kozlovsky and Lingenfelter 1979, their 
fig. 3), where E is energy per nucleon. We will use this cross-section to calculate the 7-ray 
luminosity below. 

Next, consider the shocks embedded in the stellar winds of massive stars. Owocki, 
Castor and Rybicki (1988) have shown that these shocks can have velocities (U s ) as large 
as the wind velocity (U w ) in the observer's frame. For massive O or B stars, the wind 
velocities are observed to be close to 0.01c and mass losses are of the order of 10 -5 M Q 
per yr (Bieging, Abbott and Churchwell 1989). Suppose that a fraction of i] p of the kinetic 



energy in the shock is deposited in the cosmic ray particles (rj p ~ 0.1 in, e.g., Blandford 
and Eichler 1987). The spectrum of cosmic ray particles at high energies, where E ~ pc, 
is observed to be a power law in energy E, but theoretically, the spectrum is expected to 
be a power law in momentum p (e.g., Drury 1983). We will use a spectral index a of 2.3 
in momentum in our calculation (Biermann and Cassinelli 1993). One can then write the 
kinetic energy of cosmic rays with the number density N(E)d(E) oc E~ a dE observed at 
relativistic energies , as 

V^U? = jA(J^)-*E kd (J^), (1) 

where p w is the density in the stellar wind, A (cm -3 ) denotes the abundance of cosmic 
ray protons and Ej~ is the kinetic energy. 

At energies ~ GeV per nucleon, the abundance of oxygen nuclei in cosmic rays is ob- 
served to be ~ 3.5 x 10 -3 with respect to that of the protons. In other words, the spectrum 
for oxygen nuclei, No(pc)d(pc) can be written as No(pc)d(pc) = Ao( 16 ^ eV )~ a d( 16 Q eV ), 
where Ao = 3.5 x 10 -3 A However, stellar winds of OB stars are not enriched and the 
abundance of heavy nuclei will be smaller than the standard abundance in the cosmic ray 
particles. If we we use the abundances of oxygen nuclei in the solar system (denoting it 
by Ao, s = 6.7 x 10 _4 A) then we will obtain a lower limit of the flux. We note here that 
first ionization potential effects can increase the abundance of elements of low ionization 
potential, so that Ao may be somewhat larger than Ao, s - 

The luminosity per proton (erg/s) can then be expressed as E 1 f No(pc) a(pc)vd(pc), 
where E 1 is the energy of the emitted photon (here, 4.4 MeV). The total luminosity of the 
4.4 MeV line from a region of with volume V and density p, is given by 

L 1 = n a VE 1 ( ( N (pc)a(pc) v d(pc)^ . (2) 

Here, n a denotes the density of proton in the ambient dense gas. For a typical value 
of the density, we can use the numbers that Genzel and Stutzki (1989) have tabulated for 
the physical parameters of the dense molecular gas surrounding the OB stars in Orion A 
(their Table 1). The typical density is n a ~ 10 5 cm -3 and kinetic temperature, T ~ 100 
K. 

The mass density in the stellar wind, p w , in eqn (1) can be written as p w = 
(M 2 /47rr 2 U w ) , where M is the mass loss of the star, r is the distance from the star 
and U w is the stellar wind velocity. Here the relevant volume is ~ nr 3 , where r is the 



current radius of the shock in the stellar wind. Consider a completely spherical shock 
travelling outwards through the stellar wind. The time evolution of the 7-ray line luminos- 
ity is thus an initial rise, when particle acceleration starts; a 1/t decay follows with a strong 
rise to a maximum when the shock hits the surrounding material. The radius r where we 
envisage most of the the production of 7-ray photons, that is, where the shocks meet the 
surrounding dense molecular or atomic gas, is where the ram pressure of the stellar wind 
balances the thermal pressure of the gas. In other words, (M/47rr 2 U w )U 2 ~ n a kT, which 
leads to, 

r pc ~ 0.62 n- a y 2 T^ 2 {M_^uH 2 _ 2 _, . (3) 

Here, r is in parsecs, n a = 10 5 n ai 5 cm -3 , T = 100 T 2 K, M = 10 _5 M_5 M per yr, 
and U w = 10~ 2 - 5 U w - 2 .5 c. 

Shock in the wind for OB or Wolf Rayet stars are only marginally super-Alfvenic or 
even sub- Alfvenic at the equator. The relevant Alfven speed is proportional to sin# along 
a spherical surface such as a shock in our simplified picture. The Alfvenic Mach number 
of the shock, and the density jump, are very small near the equator, and the spectrum 
is steep. Steeper spectra give rise to more 7-ray line photons and a higher ratio of 7-ray 
line to the 7r°-decay luminosity. We have calculated both luminosities as a function of the 
equatorial Alfvenic Mach number and find that for a Mach number that is below unity 
at the equator, the 7-ray line luminosity is much higher than in eqn (4), and the ratio of 
the line luminosity at 4.4 MeV to the 7r°-decay luminosity at 100 MeV is large enough to 
explain the data. In the following we use the case of the equatorial equivalent Alfvenic 
Machnumber 0.5, close to what we estimated on the basis of the nonthermal radio emission 
(Biermann, Strom, Falcke 1994); one consequence is that no particles are accelerated in 
an equatorial belt. 

Combining equations (1), (2), and (3), and allowing for the latitude variation of the 
Alfvenic Machnumber of the shock, one obtains the 7-ray luminosity at 4.4 MeV, as (using 
the cosmic abundance of oxygen nuclei as a lower limit, as discussed above) 

L 7 > 3 x 10 33 r pc n a , 5 V P ,-i (M_ 5 )(t/ s ,- 2 .5) 2 (^,-2.5)" 1 erg/s 

« 2 x 10 33 r^ 2 T" 1/2 ^ Pi _ 1 (M_ 5 ) 3/2 (£/ w ,- 2 .5)- 1/2 (t/ s ,-2. 5 ) 2 erg/s , (4) 

where r\ v = 0.1rj Pj -i. The opacity of the gas for these photons due to, for example, 
Thomson scattering, is small for reasonable parameters of the wind and the luminosity in 



the above expression is, therefore, will not be diminished. At a mean distance of ~ 0.45 
kpc for the I Orion OB association (Genzel and Stutzki 1989), and for 4.4 MeV photons, 
this means a flux of 5 x 10 -5 photons cm -2 s _1 . The observed flux in the energy range 
3 — 7 MeV is ~ 10 -4 photons cm -2 s _1 , of which about 50% is thought to be due to the 
4.4 MeV line (Bloemen et al. 1994). 

There are 56 06 through B2 stars in the I Orion OB association (Genzel and Stutzki 
1989); as they also mention, there is evidence mainly from infrared data that there may be 
a few more luminous stars with strong winds, embedded in denser gas. Therefore, the ob- 
served flux can be accommodated within the framework of the above scenario, considering 
the uncertainty in the geometry of the shocks in the wind, the actual mass loss rate and 
the shock and wind velocities. We implicitly assume that the duty cycle of shocks hitting 
the surrounding medium is fairly high; the time-scale involved is of order 10 3 yrs. It is 
worth noting that, we do not require an increased cosmic ray intensity in the clouds; we 
use only the parameters for energetic particle acceleration in stellar winds, as supported 
by radio observations of OB stars and radio supernovae (Biermann and Cassinelli 1993). 

With the cosmic abundance of oxygen nuclei in the ISM, it is easy to see that the 
contribution to the above luminosity from the collisions of energetic protons in cosmic rays 
with ambient oxygen nuclei is same as in eqn (4) . However, the width of such a line will 
be much smaller due to the small recoil velocity of the heavy nuclei in a collision with 
energetic cosmic ray protons. 

3. DISCUSSION 

Bloemen et al. (1994) noted the absence of any contribution to the luminosity in the 
3 — 7 MeV range from n° decay. It is, therefore, worth while to calculate the 7-ray 
luminosity due to decaying n° from proton-proton collisions in the above scenario. Using 
the cross-section for neutral pion production in p — p collisions given in Stephens and 
Badhwar (1981), we find that the luminosity in the pion-decay 7-ray photons at 100 MeV 
photon energy is, 

L 7)7r o «3x 10 33 r pc n a , 5 r) p _i (M- 5 )(U S - 2 . 5 ) 2 (U W -2.5)" 1 erg/s . (5) 

Bloemen et al. (1994) have discussed the emissivity of the Orion complex for energies 
more than 100 MeV, as was observed by COS-B (Bloemen et a/.1984). The luminosity of 
the cloud complex at ~ 100 MeV is ~ (4 — 7) x 10 33 erg s -1 . This is consistent with our 
derived expected fluxes. 



We have also calculated the luminosity from bremsstrahlung emitted by energetic 
electrons accelerated in the shocks. Using the detected flux in the lines as the upper limit 
of any background continuum, the electron bremsstrahlung does not, however, lead to any 
additional constraints on the ratio of electrons to protons in cosmic rays. 

It is important to consider the energy loss of the accelerated heavy nuclei in the stellar 
wind, apart from the line radiations discussed above. The process of acceleration of the 
heavy nuclei competes with energy losses due to interaction with the plasma and it is 
important that they are indeed accelerated to the required energies to be able to produce 
the observed 7-ray lines. In the denser region near the stars, the energy loss is large and 
one expects a critical length scale (the distance from the star) only beyond which the 
injection of the heavy nuclei to high energies is possible. 

To estimate the critical length scale, we consider the energy loss of the heavy nuclei 
in interactions with the thermal plasma. At the postshock temperature of T s ^ = 1.4 x 
10 7 (t/ s /10 3 kms -1 ) 2 K, the gas will be highly ionized (and therefore, ionization losses will 
not dominate). The energy loss for low energy oxygen nuclei can be written as (per 
acceleration cycle as in Biermann (1993)), 

A[3 » 9.7 x 10- 6 (t/ s ,_ 2 . 5 )- 3 (t/ w ,-2.5)" 1 M-5 r" 1 , (8) 

If we equate this loss with U s /c, the gain per cycle in first order Fermi acceleration, we 
obtain a critical length scale of r cr u ~ 10 16 (L r S) _2.5) _4 (t/ u , ) _2.5) _1 M_5 cm. In a supernova 
event the initial peak line luminosity thus scales with the fifth power of the shock velocity 
provided (a) the photons can escape, and (b) spallation does not stop the acceleration 
altogether. This luminosity arises from the shock travelling through the stellar wind. 
Thus, the nuclei can readily be accelerated as described in our scenario above. 

There is possibly a large contribution to accelerated particles from the steady stand- 
off shock of the stellar wind when it encounters the dense circumstellar material. If cosmic 
rays were injected at that shock, and cosmic rays would influence the shock structure, 
then a strong instability develops (Zank et a/.1990). In that case, the dense circumstellar 
gas is likely to be mixed into the shock region, making injection of energetic particles 
rather difficult. Therefore we emphasize that the shock travelling through the stellar wind 
and then hitting the outer boundary of the wind zone is more likely to be effective in 
accelerating particles. 

The energetic particles from these shocks ionize only a small fraction of the overall 
cloud volume, so that outside the shells a cosmic ray spectrum with a turnover, or a cutoff, 
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at low energies will be expected. The cutoff/turnover for cosmic rays in the ISM is expected 
to be ~ 30 MeV for protons from the observations of ionization rate in the diffuse clouds in 
the ISM (Nath and Biermann 1994). The ionization loss from interaction of the energetic 
particles in the wind shell limits the the grammage seen to a few g/cm 2 , otherwise the 
ionization loss would eliminate those particles which are responsible for the excitation of 
the energetic C and O atoms; beyond the shells low energy protons arising from the wind 
shocks are no longer effective. 

The production of 26 A1 in the scenario proposed may be relevant in accounting for its 
observed high abundance (Lee et al.1977; Clayton 1994); we will consider this question in 
a separate communication. 

There are various theories of cosmic ray origin that may explain the 7-ray line ob- 
servations; in particular, the works of Bykov et a/. (Bykov and Toptygin 1990; Bykov and 
Toptygin 1992; Bykov and Fleishman 1992; Bykov and Bloemen 1994), who argue for ac- 
celeration by a collection of shock waves. In the theory of cosmic ray origin proposed by 
Biermann et al. (Biermann 1993; Biermann and Cassinelli 1993; Stanev et al. 1993; Rachen 
et aZ.1993), one important argument is the acceleration of cosmic ray nuclei in SN-shocks 
that go through stellar winds. These shocks are argued to provide particle energies to 
~3x 10 9 GeV, dominated by heavy elements. The interpretation given here is thus an 
important consistency check on this proposal. 

CONCLUSIONS 

We have shown that the recently detected 7-ray lines in 3—7 MeV range by COMPTEL 
in the Orion complex can be explained in terms of emissions from energetic heavy nuclei 
accelerated in the shocks of stellar winds of massive stars as the nuclei collide with the dense 
and cool gas in the surroundings. Our calculations does not contain any free parameter 
and the luminosity is derived as a function of the ambient density, mass loss rate of massive 
stars, velocities of the stellar w inds and the shocks within. We emphasize that the stellar 
parameters used are the same as those which we need to interpret the nonthermal radio 
emission, stellar winds. 

We thus provide here a consistency check on the theory that nuclei get accelerated in 
shocks in the winds of massive stars. 
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